Quantum effects for ballistic transport in spintronic devices 
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Recent fabrication of atomic precision nanodevices for spintronics greatly boosted their perfor- 
mance and also revealed new interesting features, as oscillating magnetoresistance with number of 
atomic layers in a multilayered structure. This motivates the need to go beyond the usual theoreti- 
cal approach of semi-classical continuous layers. Here the simple tight-binding dynamics is used to 
describe quantum conduction in a multicomponent system with spin-polarized electrodes separated 
by an ultrathin and atomically coherent non-magnetic spacer (either metallic or insulating). A pos- 
sibility is indicated for obtaining a huge resonant enhancement of magnetoresistance in such device 
by a special choice of gate voltage on the spacer element. 

PACS numbers: 34.80.Pa; 73.50.-h; 73.61.-r; 75.70.Cn; 85.30.Mn 



INTRODUCTION 

In our information based society the development of 
ultra-high density storage technology is a demanding pri- 
ority. In this context, the necessity in ultra-high sensitiv- 
ity read-head devices is a great challenge from both the- 
oretical and experimental points of view. Presently the 
most promising candidates for this purpose are the mag- 
netic tunnel junctions (MTJ) made by two magnetic elec- 
trodes separated by an ultra-thin non-magnetic spacer 
and their study becomes one of the central topics in the 
fast developing field of spintronics. Since the early stud- 
ies by Tedrow and Meservey [l| on spin polarized tunnel 
conduction, an impressive progress was achieved either 
in experiment and in theory (E,@,[3| for the spin- 

tronics applications of this mechanism. The most im- 
portant recent advances are related to nano-fabrication 
of multilayered systems on atomic precision level @, 0|, 
which rises the MTJ performance up to 400% of mag- 
netoresistance and enables a breakthrough to their fun- 
damental quantum properties. It should be noted that 
the overall number of electronic degrees of freedom in 
a device like MTJ is macroscopically big which gener- 
ally suggests a quasi-classical behavior, controlled by the 
spin-dependent relaxation times or by the spin-dependent 
tunneling amplitudes. But the essentially quantum be- 
havior turns to be possible at effective separation of a 
small number (few units) of electronic degrees of free- 
dom among the macroscopically big total number as, e.g., 
the hoppings between the planes in the spacer among all 
possible hoppings in a junction, forming a partial quanti- 
zation of energy spectrum and drastically enhancing the 
sensitivity of tunnel (or ballistic) transport to external 
factors [10|| . Another natural quantization effect is the 
oscillatory behavior of conductance, either in function of 
the spacer thickness (or, more exactly, the number of 
atomic planes) and in electric field on the junction [111 ]. 



which may also allow an interesting possibility for specific 
magnetoconductance oscillations. All this needs that the 
mode mixing due to the interface roughness and intra- 
spacer defects be below the characteristic energy quan- 
tization scale, and practically requires that the spacer 
consist of few atomic planes, coherent enough. 

Consequently, the theoretical analysis of such systems 
requires a fully quantum-mechanical description, rather 

m 



than more traditional semi-classical approaches 
Up to the moment, there already exists a rather well 
elaborated theoretical base for this description, using the 
Green function formalism and rigorous ab-initio band 
calculations 
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lq as inputs to the general Kubo's 
formula. However, in many practical cases the direct use 
of corresponding algorithms leads to heavy enough nu- 
merical work, specific for each particular configuration 
and not very well suited for qualitative predictions and 
optimization of device performance. 

In this paper instead, the simple tight-binding dynam- 
ics in single-band approximation is developed, using the 
straightforward equations of motion for on-site quantum- 
mechanical amplitudes, to get a handy description of 
quantum magnetotransport in the ballistic regime (ab- 
sence of either thermal or impurity scattering) for a tri- 
layer system of spin-polarized electrodes with an ultra- 
thin and atomically coherent non-magnetic spacer. The 
motivation for our approach is an easy generalization 
to more promising device geometries (double barriers or 
double junctions, etc.) and conduction regimes (includ- 
ing finite electric field effect) which will be presented in a 
forthcoming work. This presentation is mainly limited to 
the basics of the method and to its most characteristic re- 
sults. Thus in Sec. 1 the explicit quantum wave functions 
are obtained for the 1-dimensional isolated atomic chain. 
In the following Sec. 2 the finite ID chain is inserted 
between two ID semi-infinite leads and the transmission 
and reflection coefficients for a collective electronic state 
are analytically calculated. Further, in Sec. 3 this result 
is generalized to the 3-dimensional case and the Landauer 
conductance formula [17] is used in the 3D version to 
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FIG. 1: Finite atomic chain with tight-binding amplitude t. 



yield a clear picture of basic quantum effects evolved in 
this coherent system. In Sec. 4, the important effects 
of electronic correlation are included into consideration 
using the approximation of phenomenological interface 
potentials, which foresees a more consistent treatment in 
the spirit of density functional theory. At last, in Sec. 5 
a work summary and the principal results are presented 
and commented. 



1. BASIC CHAIN MODEL 

The simplest model for transport over exact electronic 
states considers a linear chain of n identical atoms with 
single available electronic state |^) on each l-th atomic site 
and describes the single-electron dynamics in the simplest 
tight-binding approximation with (real) hopping ampli- 
tude t between nearest neighbor sites (taking the distance 
between them as unit length). 

In this coupled chain, any collective electronic state 
can be expressed as |c) = J>2i=i ci\l), with complex am- 
plitudes ci and atomic states \l) — cj|0), generated by 
the second quantization operators acting on the vacuum 
state |0). Choosing the on-site atomic energy as a refer- 
ence (e c = 0), we write the Hamiltonian operator as: 



2xu n (x) = u n+1 (x) + u n -x(x), 



(3) 



with uq(x) — 1, ui(x) — 2x. A useful trigonometric 
parametrization u/(cos#) = sin[(l + 1)9]/ sin 8 permits 
to present the general solution of Eq. [3] as: 



ui(x) 



■ [0 + l)gx] 

sin q x 



(4) 



where q x — arccosrr. Then the discrete energy spectrum 
resulting from zeros of u n [x) is explicitly given by: 



Err, = 2t cos ■ 



1 



m = 1, . 



(5) 



Now let c(x) — (ci(x), . . . , c n (x)) be the eigen- vector of 
the Hamiltonian matrix, Eq. [21 related to the eigen- 
energy e = 2tx (in what follows we mostly drop the ex- 
plicit energy arguments of amplitudes like q). Its com- 
ponents satisfy the tight-binding equations of motion 



2xci = ci+i + ci-x, for 2 < I < 



1. 



(6) 



completed by 2xc\ — C2 and 2xc n — c n -i- Since Eq. [6] 
for q/ci is just equivalent to Eq. [3] for the eigen- 

vector components can be written as: 



sm (lq x ) 

ci = — : ci. 

smq x 



(7) 



We notice that this solution also satisfies the above men- 
tioned equations of motion for c\ and c„ and provides 
the "closed" boundary conditions: 



i=i 



'1+ 



A) 



(i) 



and obtain the electronic spectrum e TO (m = 1, . . . , n) as 
the roots of the secular equation D n {e) = det(e-H^) = 

(n) 

with the corresponding Hamiltonian matrix H t J = 

(l\H^\V) =t(8i,i+i6i-i+6i,i-i6 n -i) (where 5 U > is the 
Kronecker delta and 9i — 1 if I > 0, otherwise zero). 
These determinants satisfy the recurrent relation: 



D n [e) = e£>„_i(e) - t 2 D n _ 2 {e), 



n > 2, 



(2) 



with the initial conditions Dq(e) — 1, Di{e) = e, that de- 
fine them exactly through the 2nd kind Chebyshev poly- 
nomials: D n (e) — t n u n (s/2t) Then it is convenient 
to pass to these dimensionless polynomials u n (x) as func- 
tions of the dimensionless variable x = e/2t, rewriting 
Eq. [2] as: 



co = c„+i = 0. 



(«) 



As usual, the value of c\ is fixed by the normalization 



condition, J2i \ c i( a 



1, giving finally the l-th com- 



ponent of the eigen-vector (related to the eigen-energy 

£>m 2tx r n) as: 



2 . TTml 
ci (x m ) = \l -sm 



1 



1 



Our next purpose is to consider this finite chain inserted 
into the "circuit" between two semi-infinite chain leads. 



2. TRANSMISSION THROUGH DISCRETE 
CHAIN STRUCTURE 

For a composite system of finite n-chain (in what fol- 
lows called the gate, G) between two semi-infinite chain 
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FIG. 2: Composite system of finite n-chain (gate element, G) 
inserted between two semi-infinite chain leads (source, S, and 
drain, D). The energy diagram shows the on-site energy levels 
(dashed) for i-th element (i = s,g,d) and the Fermi level 
(dot-dashed) whose crossings with the continuous S- and D- 
dispersion curves define the wave numbers for incoming (q s ), 
reflected (— q 3 ) and transmitted (qd) parts of the Fermi state. 
Notice that the Fermi level generally does not match any of 
the discrete levels (solid) in the central (G) element. 



leads, S (source) and D (drain) (Fig. [2]), the tight- 
binding Hamiltonian, Eq. [JJ is extended to: H = 
h s + h d + h g + h s 9 + hP d where: 



oo 

h S = ^ [ £ ssl*l + *« (4^+1 + n - C - 



1 = 1 



h 3 
h sg 



n 

[ e a9i9l + t g (gjm+i + h.c.) 



/=1 



t sg (s\gi + h.c. 



(9) 



including the respective on-site energies e, (i = s,d,g) 
and hopping parameters t, (i = s,sg,g,gd,d), while the 
operators h d , h ad are analogous to h s , h sg with the formal 
change of indices s — > d. For this macroscopic system, 
the energy spectrum includes continuous S- and D-bands 
£i iQ = e i + 2tiCosq, i = s,d and possibly discrete G- 
levels outside these bands. The collective eigen-state for 
a given energy e can be found from the equations of mo- 
tion that generalize Eq. [6l We denote s/, <?/ or di the re- 
spective local amplitudes of the wave function and define 
the dimensionless dynamical variables Xi = (e — ei) /2ij 
(i = s, <?, d). Let the S-amplitude be a sum of an incident 
wave of intensity 1 with the wave number q s — arccos x s 
and a reflected wave with certain amplitude R and the 
wave number — q s : 
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(10) 



(for regressive order of sites I in S), and the D-amplitude 
present a transmitted wave with certain amplitude T and 
the wave number qd = arccos Xd- 



(11) 



Eqs. [10] and [TTJ refer to one of fundamental solutions 
for given e (besides that where the incident and reflected 
waves belong to D and the transmitted one does to S). 
These forms automatically satisfy the equations of mo- 
tion within S and D: 



2x s si = si-i + sj+i, 2x d di = + di+i 



(12) 



(for I > 2) , while the pairs of equations on the S/G and 
G/D interfaces: 



2si cos q s 
uigi 



and 



2d\ cos qd — 



u\g„ 



. ^sq 

S2 + -r-gi, 

92 + -r-Si, 



«2 + — 9n, 
td 



(13) 



g n -i 



to 



(14) 



are the discrete analogs of usual boundary conditions for 
continuous wave function and its derivative jifj. They 
permit to express the terminal pairs of G-amplitudes 
through the asymptotic parameters R, T, q s and qd- 



.91 



t sg 

td rp 

x 1 ' 

tgd 



(1 + R) , .92 



' sg 



[ui - 7* + (ui - 7 S ) R] 



-x = -^{u x - ld )T, (15) 

tgd 



with the interface parameters 7 S = e lq3 t sg 2 /t g t s and 
7d = e zqd t g d 2 /tgtd- The polynomials ui = ui(x g ) are 
formally the same as given by Eq. [4] with the energy 
argument x g — (e — e g )/2t g . But the energies e of our 
main interest for the transport processes are those close 
to the Fermi energy e-p which is generally not an eigen- 
value, Eq. for the isolated G-element. Therefore the 
transient "momentum" q g — arccos x g (not necessarily 
real) breaks down the closed boundary conditions, Eq. 
for G and thus enables continuity of quantum states 
along the composite system. Next, using Eq. [6] for this 
element in the form: 
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u%9l =91+1+91-1, (16) 
it is possible to interrelate the terminal G-amplitudes: 



g n ~l = U n -29l - -J^U n -3Sl 



9n = Un-l9l 



t, 



-U n - 2 Sl. 



(17) 



Then, Eqs. [15] and [TTl yield two independent relations be- 
tween the coefficients R and T. Those are readily solved 
to give: 



R{x S , Xgj Xd) . 



T (X S , Xg , X(l) 



(18) 



where the resonance properties result from the denomi- 
nator: 
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FIG. 3: Transmission coefficient \T\ in function of the on-site 
energy e g in the gate element of the composite chain system, 
for the choice of its parameters e s = —0.4 eV, Ed = —0.8 eV 
(relative to the Fermi energy), t s = td = 0.5 eV, t g = t sg — 
t g d = 0.25 eV and n = 5. The shadowed areas indicate the 
(Stoner shifted) continuous bands, S (light grey) and D (dark 
grey). 



Dn [JEs 7 %g s Xd) 



{is + 7s) u n-l + lsldUn-2, (19) 



with the relevant variables Xi as arguments of complex 
factors 7i and real polynomials ui, and D n {x s ,x g ,Xd) = 
D n {x s + ir,x g ,Xd). Since, in the considered ID case, all 
Xi(s) = (e — Si) /2ti are defined by the single energy vari- 
able s, the coefficients R and T can be also defined as 
functions of energy: R(e) = R(x s (e) 7 x g (e),Xd{s)) and 
T(e) = T (x s (e),x g (e),Xd(s)). It is important to notice 
that the result of Eqs. 1 18) 191 is just analogous to that ob- 
tained with the Green function techniques 0] , the factors 
7 S and 7^ playing the role of interface Green functions of 
Ref. 0- A typical behavior of the transmission coeffi- 
cient |T(e)| 2 is presented in Fig. [31 It shows n transmis- 
sion resonances generated by n discrete energy levels of 
the G-element (by n atoms in the chain) as they are pass- 
ing over the Fermi level within the mutually displaced en- 
ergy bands. The displacement can be due, for instance, 
to the Stoner splitting between majority and minority 
subbands of oppositely polarized S- and D-elements (see 
also Sec. ). Notice that the resonances become sharper 
as the levels approach the band edges, and the maximum 
transmission in the asymmetric S-D band configuration 
is not limited to unity. This coefficient enters the Lan- 
dauer formula [TtJ for the ballistic conductance through 
the ID composite system (in zero temperature limit): 



G=-\T(e F )\ 2 , 



(20) 



with the Fermi energy ep. Now, allowing the S and 
D chains to support spin polarized subbands £i,q >a — 



Ei.a + 2tiCOsq (where Ei >a = E{ — crAi, a = ± are 
the majority and minority spin indices and A, are the 
Stoner splitting parameters for i = s,d), we can intro- 
duce the energy and spin-dependent variables Xi l(T (s) — 
(e — Ei t(T ) I (2ti), i = s,d, for in- and out-channels and ob- 
tain from Eq. [20] the spin-dependent conductance values 
G a , a t = (e 2 /h)\T(x s ^(E F ),x d ^,(£ F ),x g (£ F ))\ 2 . Fi- 
nally, the (maximum) magnetoresistance is defined as 
usually through the difference between the conductance 
values Gp = G+,+ + G_._ for parallel and Gap — 
G + ,_ + G + ,_ for antiparallel S/D polarization: 



MR=^-1. 
Gap 



(21) 



Although the state-of-the-art technology alr eady permits 
development of such genuinely ID devices [19( and the 
resonance behavior like that in Fig. [3j (different from 
the known quantized conductance steps vs voltage bias) 
can be directly sought in them, it is of major practical 
importance to generalize the above treatment for a more 
realistic multilayered structure and this will be done in 
the next section. 



3. 3-DIMENSIONAL MULTILAYERED 
STRUCTURE 

Passing from ID composite chain to multilayered (and 
spin polarized) 3D lattice structure as shown in Fig. 01 
we extend the indexing of site operators from §1 , di and gi 
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FIG. 4: Real multilayered structure where the current / flows 
through two ferromagnetic electrodes, S and D, separated by 
a non-magnetic spacer G and its model by the composite 3D 
system where a finite n-plane spacer is inserted between two 
semi-infinite leads. 

to Si,m,a, ^,m,cr and <?i, m ,cr, where m runs over N sites in 
the lih plane and a is ±. Our strategy in this case relies 
on the conservation of the transversal quasi-momentum 
k = (k x , k y ) in the transitions between the planes [IRIilj . 
From the experimental point of view, this requires perfect 
interfaces that are only reachable with advanced molec- 
ular beam epitaxy (MBE) techniques [12] • To describe 
the situation where k is a good quantum number for inde- 
pendent ID-like longitudinal transport channels, we pass 
from the site operators to the "planar wave" operators. 
Thus, for the Ith plane in the S element, the latter oper- 
ators are defined as: 



turn subbands £i,k,cr == £i,a- + %U (cos k x + cos k y ) ,i = s,d 
and £(,,k,er = £ g + 2i s (cosk x + cosk y ). The equations of 
longitudinal motion in terms of the "planar wave" am- 
plitudes s/,k,cr, di.\L,a and gi^a (for given energy e of the 
collective state) are obtained in analogy with the ID case. 
Thus, in the leads S and D (beyond the interfaces, at 
I > 1), they are analogs to Eq. [321 

2a; s ,k,o-Si,k,CT = Si-l,k,<r + S2 + i,k,cr, 

2aJd,k<^,k,cr = d;-i,k,<T + rf/+i : k : <r, (23) 

with Xi,k,<r = {s — £i,k,a) I (2ti), i — s.d, while in the 
spacer G (at 1 < I < n), we have in analogy with Eqs. [7| 
and [TBI 

2a; g ,k3i,k : <j = .gz-i.k.o- + g;+i,k l( r, (24) 

with Xg k = (e — £ ff .k) / (2i ff ). Also the equations for in- 
terface amplitudes: 



2a; s ,k,(TSk,l,o' — Sk : 2,cr + — Sk.l.CT; 

o i ^ s 9 

£%g,\i,crQ\t,\,a — 5k,2 : cr + — Sk,l,<r, 

t g 

^■^g.k.,cr9k.,n.cr — ^k,n — 1,ct "T" ~~ ^k,l,CTi 

t a 

2Xd,k,o-rfk,l,CT = rfk,2.o- + -J—gk.l.cr, (25) 

td 

are analogous to Eqs. [HI [14] and [HI The next deriva- 
tion, in full similarity with the ID case, leads to the full 
dispersion laws in the leads £i^ q ,<j — £i,k,<j + %U cos q 
(for i = s,d) and to the final conductance formula for 
particular in-out spin channels: 



J ; , k , CT - ~y=Y e ikm Si, m , CT , (22) 

v m 

and, similarly, the "planar wave" operators di^ a 
and f^k.er for D and G elements are written. The 
related extension of the Hamiltonian is: H = 
Ek, CT (k„ + H,a + K* + K% + Ka) > where the par- 
ticular terms are analogous to those in Eq. [9] with the 
change of all the site operators by the "planar wave" ones 
and all the on-site energies £j by the transversal momen- 



G ™' = \ 12 l T -'(^F,k)| 2 . (26) 

Here the transmission coefficient depends on the rele- 
vant variables a 7 a',e and k accordingly to: T acr '(e, k) = 
r(g Sj k, CT , q g ,k, qdM,tr') with q l ^. a = arccosa^ka for i = 
s,d and q g ^ — arccosa^k, and the sum in k is re- 
stricted to the "permitted" range K, such that simul- 
taneous equalities e s ,k,g„ = £d,k,q d = £f are possible for 
certain real q s and qd- In more detail, the latter condition 
is expressed as: 



max {—2, max [x Sl(r (ef) , ^d.a 1 (£f)] — 1} < cos k x + cos k y < min {2, min [x S)Cr (ep) , Xd.o 1 (£f)] + 1} , (27) 

I 

fully defining the integration procedure (in the limit of continuous k). Then, seeking for optimum performance 
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of the model MR device from Eq. [211 it is of interest 
to evaluate it as a function of the system parameters, 
mainly the number of atomic layers in the gate n and 
the on-site energy level of the gate e g (which can be pos- 
sibly controlled through the gate bias). Also, variation 
of the latter parameter from positive to negative values 
permits to model in our approach the passage from the 
tunnel magnetoresistance (TMR) to giant magnetoresis- 
tance (GMR) regime in a unified way. 

The following numerical work can be oriented accord- 
ingly to some evident qualitative arguments. The varia- 
tion of the integrand in Eq. [26] is mainly controlled by 
that of the polynomials ui (:%,k) in the denominator of 
Eq. [T9l As seen from the explicit Eq. [4] they are oscil- 
lating if |x S) k| < 1 (that is, the sampling point e 9i k in the 
G-spectrum close enough to the Fermi energy ep) and ex- 
ponentially growing if \x g \ > 1 (remote enough e g ^ from 
£f)- Therefore, the conductance is generally expected to 
oscillate (either in e g and in n) as far as the level e g is 
close enough to e-p (which can be compared to the GMR 
regime) and to exponentially decay at e g far enough 
from £f (a generalized TMR regime). The latter decay 
should asymptotically tend to MR(n) oc exp (— nx m i n ) 
with Xmin = min ke j<- |i%k| at n > 1. 

In the latter case, the direct calculation by Eqs. [21] 
may result in Gp and Gap both exponentially small 
but the latter yet much smaller and thus in (arbitrarily) 
huge MR values. However, they should not be physically 
attainable, taking into account that the real multiband 
electronic structure of transition metals always includes 
some additional conduction channels, for instance due 
to the s-bands, whose tunnel contribution slower decays 
than that of d-bands and is almost spin independent. 
Therefore it should dominate the transport in the indi- 
cated regime and make the real MR exponentially small. 
A simple phenomenological account of this mechanism 
in the considered single-band model can be done by in- 
troducing a certain spin-independent term Go into either 
Gp and Gap values: 



G P = G+++G—+G , G A p = G + -+G-++G a , (28) 
to present the MR formula, Eq. [21] as 



MR = 



G_i 



G__ -G+_ -G_ 



G- 



G 



(29) 
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FIG. 5: Magnetoresistance of a FM/NM/FM junction with 
the model parameters: e a ,+ = £d,+ = 1-32 eV, e s ,- = £d,- = 
3.36 eV, t 3 = t d = tg = -0.6 eV (like those from Ref. [H])' and 
7s, d = 0.5 in function of the number n of spacer layers at fixed 
values of e g . Notice the exponential decay in the TMR- like 
regime either at the highest e g = 4 eV and the lowest e 9 = — 1 
eV and a strong enhancement with emergence of oscillatory 
behavior at intermediate e g ("shallow band" regime). 
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FIG. 6: MR vs the on-site energy s g of the gate element for 
a junction with the same parameters as in Fig. [5]and various 
numbers of atomic planes in the gate element, n = 2, n = 4 
and n — 5. Compare the resonance peaks in the "shallow 
band" regime with those in the ID case of Fig. [3] 



It is just this formula that is used below for all practical 
MR calculations. Thus, using the band structure pa- 
rameters: £ s .+ = Sd,+ = 1-32 eV, £ S! _ = £d,- = 3.36 
eV, t s = td = t g = —0.6 eV and j Si d — 0.5 (a rea- 
sonable single-band fit for the real Fe band structure, 
see 0,0, El) and choosing for simplicity the constant 
value Go = 0.1e 2 /h, we find that the MR behavior vs 
n indeed changes qualitatively at different choices of e g 



(Fig. [5]). The TMR- like behavior with fast exponential 
decay appears either at high enough gate level, e g > 6t g 
(which can be compared to a "positive" barrier in the 
continuum approximation), and at low enough e g < — 2t g 
(a "negative" or "hole" barrier), but it develops GMR- 
like oscillations with greatly increasing overall MR ampli- 
tude at the intermediate e g values (which can be called 
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the "shallow band" regime). The oscillating behavior 
is in a qualitative agreement with that experimentally 
observed for MgO moderate tunnel barriers between Fe 
electrodes [§], except for stronger first oscillations than 
in the data. However, it will be shown below that these 
strong oscillations are effectively moderated when taking 
into account the specific interface effects between metal 
and insulator layers. The most notable feature of the 
calculated MR is its amazingly high maximum value, of 
the order of 3000 %, indicating a huge potentiality of the 
quantum coherent conduction regime. 

For the same choice of parameters, the calculated de- 
pendencies of MR vs e g (at fixed values of n = 2, 4, 5) 
are shown in Fig|6l In concordance with the above con- 
sidered MR(n) behavior, they practically vanish beyond 
the range of intermediate e g and display a finite number 
of resonance peaks within this range (reminiscent of n 
ID resonances in Fig. [3]), reaching the same highest or- 
der of magnitude in the "shallow band" regime. These 
very high values in the present tight-binding approach 
contrast with the known result for the model of almost 
free electrons on the continuous rectangular barrier [io| . 
where MR reaches zero minimum at low barrier height. 
As yet, the MR(e g ) dependence was only studied exper- 
imentally for Al-0 spacers (2^|, possessing most probably 
polycrystalline or amorphous structure and high enough 
e g , so it could be of interest to try it also with epitaxial 
MgO spacers and possibly with those spacer materials 
that can realize the "shallow band" condition. 



4. INTERFACIAL EFFECTS 

In this section, we will discuss the interfacial effects 
present at the metal/insulator or metal/non-magnetic- 
spacer interfaces. This is motivated by the anal- 
ogy with the well known case of Schottky barrier at 
metal/semiconductor interfaces, leading to such interest- 
ing physical effects as band bending [27||. It is known 
from X-ray and ultraviolet photoemission spectroscopy 
(XPS and UPS) that some charge transfer effects also 
appear at the metal-insulator interface, leading to for- 
mation of an interfacial charge-dipole whose mag nitude 
is defined by the localized states at interfaces [28|. Since 
this dipole directly affects the efficiency of tunneling, it 
is also important to evaluate its effect in the magnetore- 
sistance. 

The best treatment of this problem is to introduce self- 
consistently a charging energy (8, commonly called the 
built-up voltage) due to a charge accumulation in the 
framework of the density functional theory (DFT). This 
is going to be done in the future work, but at the moment 
we will develop simple analytic formulas to take into ac- 
count these interfacial effects qualitatively. Despite of 
its simplicity, the model can exemplify in which way the 
formation of charge dipoles affects the magnetoresistance 
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FIG. 7: Schematic representation of the interface charge- 
energy, 5, created by a charge accumulation in the S/G and 
G/D interfaces, as a simplified description of the true self- 
consistent behavior. 



ratio. 

We go on using the same model of Sec. but consid- 
ering extra charge energies ±8 on the sites pertaining to 
the two atomic planes on both sides of each interface (see 
Fig. as an approximation for more realistic charge and 
potential distributions around interfaces, obtained by nu- 
merical DFT calculations The (5-perturbation re- 
sults in new boundary conditions and, as a consequence, 
in a new transmission coefficient. In this approximation, 
there is no changes in equations of motion within the 
particular elements (S, D and G), but new pairs of equa- 
tions appear at the S/G and G/D interfaces, involving 
the charge energy 8: 

(2cosq s + 6/t s ) si = s 2 + (t sg /t s )gi, 

(Xg - 5/tg) .91 = g 2 + {t S g/tg) SX, 

(2 cos q;, + 5/td)di = d 2 + {tgd/td) g n , 

(x g - 5/tg) g n = + (tgd/tg) d\. (30) 

These boundary conditions allow to re-calculate two 
terminal G-amplitudes in function of the parameters 
R, T, q s ,qd- Interconnecting these terminal amplitudes 
leads to the transmission formula, Eq. [TBI but with the 
modified denominator D n $ = A n — B n + C n where: 




' ( U n -1 ~ ~ U n-2 | ■ 



C n = lsldUn-2- (31) 

It is easy to see that D ni s — > D n in the limit of 5 — > 0. 
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FIG. 8: Effect of the interface charge energy on the MR(n) 
dependence. The system corresponds to the parameters as in 
Figs. 02 [6] with e g = 1 eV and S varying from to 0.4 eV. 



wells. This idea of charge energy induced resonances is 
corroborated by the calculated sharpening of peaks just 
in the AP conductance. Amazingly high peak MR val- 
ues, reaching ~ 3000% for a reasonable choice of 5 ~ 0.4 
eV (similar to the numerical estimate for Fe-MgO inter- 
faces, [l4|). should motivate fabrication of new devices 
with the choice of such spacer materials as semiconduct- 
ing (Ge, Si) or semimetalic (Sb, As). Though the peak 
value may be obviously decreased under the effects of 
electron-electron, electron-phonon, and electron-magnon 
interactions, finite temperature and disorder, it expresses 
one of the principal results of this work, demonstrating 
that the highest possible MR value should be reached in 
the shallow-band regime for non-magnetic spacer by ad- 
justing to the strongest resonance condition, once elec- 
tronic coherence is assured. 



5. CONCLUSIONS 
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FIG. 9: Effect of the interface charge energy on the MR(e a ) 
dependence. The system corresponds to the parameters as in 
Figs. [HE] with n = 5 and 8 varying from to 0.4 eV. 

The MR defined from Eqs. 1291311 in function of the 
number n of gate atomic planes and of the gate voltage 
e g for tree values of the interface potential 6 are pre- 
sented in Figs. [H [3 The obtained softening of first os- 
cillations makes these curves more similar to the exper- 
imental observations [3|. An unexpected result is that 
the effect of an extra barrier due to the charge energy 
can yet reinforce the calculated MR peak in the shallow 
band regime, though reducing the values in the TMR 
regime at higher barrier height. Obviously, the charge 
energy barrier reduces the conductance (either in P and 
AP configurations), but the MR enhancement is mainly 
due to a much stronger reduction of the AP conductance. 
Apparently it results from the wave function localization 
caused by coherent resonances in the interfacial potential 



In this work a theoretical approach was developed 
to fully coherent spin-dependent quantum transport in 
nanolayered magnetic junctions, using single-band tight- 
binding model with explicit equations of motion for wave- 
function amplitudes. The analytic solutions for the trans- 
mission and reflection coefficients were generalized for 
a 3-dimensional magnetic junction structure. The sim- 
ple zero temperature calculations have revealed the most 
pronounced enhancement of the magnetoresistance in 
the "shallow band" regime, defined by low gate voltages 
(Fig. [6]) . Another important feature for this gate voltage 
regime is the calculated oscillatory behavior of MR with 
the number of atomic planes in the spacer. In support of 
these theoretical predictions, the calculations also reveal 
that the oscillatory regime starts already at moderate 
gate voltages, e g ~ 2 eV. This agrees rather well with the 
experimental observation by Yuasa [1] of clear MR oscil- 
lations at low enough gate voltage barrier e g ~ 0.4 eV 
in a Co/MgO/Co structure. So it is concluded that the 
best MR values for a quantum magnetic junction could 
be reached using shallow band materials for spacer layers, 
the pos sible candidates sought between transition metals 
(Cr [3l| in junctions of the type Fe/Cr/Fe or Zn in junc- 
tions of the type Co/Zn/Co), semiconductors (Ge, Si), 
or semimetals (Sb, As). Finally, the important effect of 
charge build-up (Sec. ) on the junction interface was also 
considered and, though in a simple phenomenological ap- 
proach, a possibility is shown for even stronger enhance- 
ment of the magnetoresistance in presence of the extra 
charge-energy for the same shallow band regime, empha- 
sizing again the promise of using the low e ff -materials. 
To verify these model predictions, the future work must 
include various realistic effects, such as those of finite 
temperature and self-consistent density functional theory 
(DFT), to adopt also the multiband electronic structure, 
spin-transfer processes, and non-linear conduction. 
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